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Abstract 

We study phase separation in two dimensions in the scaling limit below criticality. The 
general form of the magnetization profile as the volume goes to infinity is determined exactly 
within the field theoretical framework which explicitly takes into account the topological 
nature of the elementary excitations. The result known for the Ising model from its lattice 
solution is recovered as a particular case. In the asymptotic infrared limit the interface 
behaves as a simple curve characterized by a gaussian passage probability density. The 
leading deviation, due to branching, from this behavior is also derived and its coefficient 
is determined for the Potts model. As a byproduct, for random percolation wc obtain the 
asymptotic density profile of a spanning cluster conditioned to touch only the left half of the 
boundary. 



1 Introduction 



Boundary conditions notoriously play an important role in the theory of phase transitions. For 
a system of ferromagnetic spins taking discrete values, a pure phase of type a with translation 
invariant spontaneous magnetization below the critical temperature T c can be selected fixing all 
boundary spins to the value a and then sending the boundary to infinity. On the other hand, 
if the spins are fixed to a value a on the left half of the boundary and to a different value b on 
the right half, a pattern of phase separation between phases of type a and type b is expected 
in the large volume limit below T c , at least away from an interfacial region anchored to the 
points of the boundary where boundary conditions change from a to b. The properties of the 
phase separation and the notion of interface have been extensively studied both in two and three 
dimensions through rigorous [I], exact [2] and approximate [3] methods. The most advanced 
analytic results are available in two dimensions, where the exact asymptotic magnetization profile 
has been obtained for the Ising model [H [5] exploiting its lattice solvability at any temperature. 
This result shows in particular that in two dimensions the Ising interface has middle point 
fluctuations which diverge as the square root of the volume, a property previously proved for 
low temperatures in [6]. The result for the Ising magnetization profile also admits a simple 
interpretation in terms of passage probability of the interface through a point [7J. No exact 
result for the magnetization profile is available in three dimensions. 

In this paper we use field theory as a general framework for the study of phase separation 
in the scaling limit below T c , for any two-dimensional model possessing a discrete set of ordered 
phases and undergoing a continuous phase transition. For a strip of width R we derive the large 
R asymptotics for the magnetization profiles along the longitudinal axis in the middle of the strip 
and show that a generalization of the Ising result holds whenever the surface tension between the 
phases a and b cannot be decomposed into the sum of smaller surface tensions. The formalism 
explicitly illustrates the role played by the topological nature of the elementary excitations 
(domain walls), which for a discrete set of ground states is peculiar of the two-dimensional case. 
The interpretation in terms of passage probability holds in general, with subsequent terms in 
the large R expansion accounting for the emergence of an interfacial region with finite width in 
a way that can be understood through renormalization group considerations. 

The trajectories on the plane of the domain wall excitations of the field theory are naturally 
interpreted as the continuum limit of the boundaries of clusters made of nearest neighbors with 
the same value of the spin. In the last years the scaling properties of cluster boundaries have 
been extensively studied at criticality in the framework of Schramm-Loewner evolution (SLE, 
see e.g. [5] for a review); the application of SLE methods to the off-critical case, on the other 
hand, is up to now much more limited (see [3 [10]). The renormalization group interpretation of 
our results below T c is that the cluster boundary connecting the two boundary changing points 
on the edges and the interfacial curve between the two phases coincide as R — > oo, and then 
have the same gaussian passage probability density; when R decreases, the interfacial region 
with finite width emerges via branching of the interface and formation of intermediate clusters. 
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We write down the leading term associated to branching and determine its coefficient for the 
(/-state Potts model. When the surface tension between the phases a and b is decomposable 
the formalism leads to a multiple interface description, the Ashkin- Teller model providing an 
interesting example of this type. 

The paper is organized as follows. In the next section we introduce the field theoretical 
formalism and derive the large R results for the magnetization profiles. Section 3 is then devoted 
to the interpretation of the results and to the discussion of the interface structure. The specific 
cases of (/-state Potts and Ashkin- Teller models as well as an application to percolation are 
finally discussed in section 4. 

2 Field theoretical results 

Consider a ferromagnetic spin model of two-dimensional classical statistical mechanics in which 
each spin can take n discrete values that we label by an integer a = 1, 2, . . . , n. To be definite 
we refer to the case in which the energy of the system is invariant under global transformations 
of the spins according to a symmetry group; the spontaneous breaking of the symmetry below 
a critical temperature T c is responsible for the presence on the infinite plane of n translation 
invariant pure phases; the phase of type a can be selected starting with the system on a finite 
domain with boundary spins fixed to the value a, and then removing the boundary to infinity. 
We denote by Z a and (• • • ) a the partition function and the statistical averages, respectively, in 
the phase a. 

We consider the scaling limit below T c , described by a Euclidean field theory defined on the 
plane with coordinates (x, y); it corresponds to the analytic continuation to imaginary time of a 
(l+l)-dimensional relativistic field theory with space coordinate x and time coordinate t = iy. 
This theory possesses degenerate vacua |f2 a ) associated to the pure phases of the system. In 
1+1 dimensions the elementary excitations will correspond to stable kink states \K a b{9)) which 
interpolate between different vacuaj] |fi a ) and |fi{,); the rapidity 9 parameterizes the energy and 
momentum of the kinks as (e, p) = (m a b cosh 8, m a b sinh 0), m a b being the kink mass. In general, 
connecting |f2 a ) and |fi&) requires a multi-kink state \K aai (6i)K aia2 (92) ■ ■ ■ K an _ 1 b(6 n )) passing 
through other vacua; we call adjacent vacua two vacua that can be connected through a single- 
kink excitation. There can be kinks with different masses connecting two adjacent vacua \Q a ) 
and | fife); in such a case the notations \K a b{9)) and m a b refer to the kink with the lowest mass, 
which is leading in the large distance limits we will consider. 

Consider now the scaling limit on a rectangle (Fig. [1]) with horizontal sides of length L and 
vertical sides of length R (the origin of the coordinates is taken in the center of the rectangle), 
with the following choice of boundary conditions (boundary conditions of type ab): the boundary 
spins are fixed to a value a for x < 0, and to a different value b for x > 0. Let us denote by Z a b 
and {■ ■ - )ab the partition function and the statistical averages for the system with this choice of 

1 In a kink state \K a b(0)) the spin field a(x, y) takes the value (a) a as x — > — oo and the value (<r)i, as x — > +oo. 
Propagation of the kink in Euclidean time y = it corresponds to domain wall configurations for the lattice model. 
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Figure 1: The scaling limit below T c of a ferromagnet is considered on a rectangle with boundary 
spins fixed to take the value a for x < and a different value b for x > 0. The case L — > oo and 
-R S> £ is considered throughout the paper. 



boundary conditions. We consider the limit L — > oo and want to study properties of the system 
as a function of the width R of the resulting infinite strip, focusing on the asymptotic limit in 
which R is much larger than the correlation length £. 

Within the field theoretical formalism the boundary condition at time t switching from a to 
b at a point xq is realized by a boundary state that we denote by \B ab (xo; t)). This state can be 
expanded over the basis of asymptotic particle states of the relativistic theory. The change of 
boundary conditions at the point xq requires that kink excitations interpolating between |f2 a ) 
and | rifo) are emitted/absorbed at that point. Then, if \Q a ) and \Q b ) are adjacent vacua, the 
boundary state has the form 



\B ab (x ;t)) = e 



-itH+ixoP 



df) 

^f ab (e)\K ab (d)) + ■ ■ ■ 
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where H and P are the energy and momentum operators of the (l+l)-dimensional theory, f a b{&) 
is the amplitudj§ for the kink to be emitted at the boundary changing point, and the dots 
correspond to states with total mass larger than m a b- The partition function we are considering 
can be written as 

Z ab (R) = (B ab (x ;iR/2)\B ab (x ;-iR/2)) = (B ab (0; 0)\e~ RH \B ab (0; 0)) (2) 
and, as a consequence of ([JJ, has the large R asymptotics 



—m a hR. cosh 9 \ i „-m ab R 



Zab(R) ~ / 7T e~ m ^ co ^\fab(0)V ~ -7==i . (3) 

The specific interfacial free energy, or surface tension, is given by 

v V 1 1 Z ab{ R ) f A\ 



2 We use the normalization (K a b(0)\K a , v {6')) = 2n5 aa ,5 w 5(6 - 6') 
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where Z a (R) is the partition function for uniform boundary conditions of type a on the strip. 
Since the lowest mass state entering the expansion of the boundary state \B a (t)) for uniform 
boundary condition is the vacuum |O a ), Z a (R) tends to (fi a |fi a ) = 1 as R — > oo, so that 

Safe = m ab . (5) 

If \Q a ) and are not adjacent vacua the expansion of the boundary state \B ab (xo\t)) starts 
with a multi-kink state, and the corresponding surface tension is a sum of surface tensions 
between adjacent vacua; we defer to section 4 the illustration of this case. 

Let us denote by ex a generic component of the spin field, omitting for the time being the 
index which in general labels the different components. The magnetization profile along the 
horizontal axis in the middle of the strip for ab boundary conditions and a and b adjacent 
phases is 

(a(x,0)) ab = ^( J B a6 (0;0)|e-^/ 2+ " p ( r(0,0)e-^/ 2 - ia;P | J B Qb (0;0)) (6) 

„ _L y |^/*(0)/(^)(J^(0)|ff(O,O)|^ 

the last line being the large R limit obtained from (pQ). The matrix element of the spin field it 
contains is related by the crossing relational 

(K ab (0)\a(O, 0)\K ab (9')) = F° ba {9 -9' + iw) + 2ir5(9 - 9')(a) a , (7) 

to the form factor 

KbaiPl ~ fc) = {^aW{m\Kab{ei)KUe 2 )) , (8) 

where the vacuum expectation value (a) a appearing in the disconnected part is the spontaneous 
magnetization in the phase a on the infinite plane. When 9i — 02 = in the kink and the anti-kink 
in ([S]) have opposite energy and momentum and can annihilate each other. In 1+1 dimensions 
these annihilation configurations produce in general simple poles that have been characterized 
for general fe-particle form factors in integrable field theories (see in particular [H]). For k = 2, 
however, integrability plays no role in the determination of the residue, which for the case of 
kink excitations reads [12] 

Rese =ln F^ ba (e) = i [(a) a -(a) b }. (9) 

For R — > oo the integral in ([6]) is dominated by small rapidities and the leading contribution 
can be written as 

? r°° i 

(<T(x,0)) ab ~(*) a + —[(a) a -(a) b ] J de _— e -m ah m^+im abX 9^^ (1Q) 

where we used ([3]), ((7|) and Q, 6L = 9 — 6', and we integrated over 9 + = 9 + 9'. The last 
integral is regularized moving the pole slightly above the real axis, so that the usual relation 



3 In field theory 'crossing' a particle from the initial to the final state (or vice versa) involves reversing the sign 
of its energy and momentum. Within the parameterization introduced above this amounts to a in rapidity shift. 
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(x - iO) 



-i 



iir5(x) + p.v. x 1 finally gives 




(11) 



where the principal value of the integral in (jlOp has been expressed in terms of the error function 
erf(x) = Jq drj e - ^ 2 (see e.g. [IS])- The same result can be obtained differentiating (fTTJj) with 
respect to x in order to get rid of the pole, and then integrating the result of the integral over 
6L with the condition (er(+oo, 0)) a b = {&)b- 

For the Ising model ({cr) a = — the result (jlip coincides with the scaling limit of that 
obtained from the lattice in [H \5\ . Even in our more general setting, it shows that for R — > oo 
(a(a(m a bR)^ /m a b,0)) a b tends to the pure values (a) a or (a)b for j3 > 1/2 and a negative or 
positive, respectively, and to the average value ((a) a + (a)b)/2 for (3 < 1/2. 

The result (|lip is produced by the leading term in the small rapidity expansion of the 
emission amplitude in ([T]) and of the matrix element ([7]) . More generally, for the latter we write 



with c ab given by Q. As for the emission amplitude, it satisfies f a b{6) = fba(—&) as a 
consequence of reflection symmetry about the vertical axis. In any model in which a and b play 
a symmetric role we will have f a b{8) = fba{@) and f a b(&) = /afe(O) + 0(6 2 ). Then it is easy to 
check that the next contribution to {a(x,0)) a b produced by the small rapidity expansion is 



Notice that the error function in (jlip is leading with respect to (|13j) as R — > oo for x ~ 
(m a bR)P /m a b with f3 > 0; the two terms are of the same order for (3 = 0. 

3 Passage probability and interface structure 

The results of the previous section allow an interpretation based on renormalization group and 
probabilistic considerations. It is clear that the problem has two length scales: the correlation 
length £, proportional to the inverse of the kink masses, which is the scale of the fluctuations 
within the pure phases, and the width R of the strip, which sets the scale at which we observe 
the system with ab boundary conditions. In an expansion around R/£ = oo the leading term 
corresponds to the crudest description of phase separation in which all short distance features are 
washed out and one is left with two pure phases sharply separated by a simple curve connecting 
the two boundary changing points (Fig. [2^,). Hence, the notion of curvilinear interface naturally 
arises in this limit and can be formulated directly in the continuum. It is clear, however, that this 
picture cannot hold true for finite values of -R/£, and that one needs to switch from the notion 
of sharp separation and curvilinear interface to that of an interfacial region (or thick interface) 
with a width which diverges as the correlation length when the critical point is approached: such 




(12) 




(13) 
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(a) 



(b) 



(c) 



(d) 



(e) 



Figure 2: Some configurations of the interfacial region. Lines correspond to domain walls 
between different phases. 



a divergence simply reflects the fact that phase separation disappears at criticality. Within the 
large R/£ expansion the leading deviations from the simple curvilinear picture are expected from 
effects such as branching and recombination as well as self-intersection of the curved (Fig. [2}>d). 
These effects appear in the expansion through insertions (delta functions) localized on the curve 
separating the two pure phases. 

According to this discussion the large R/£ expansion for the magnetization at a point x on 
the axis y = when the interface passes through a point u on this axis can be expected to start 
as 

a ab (x\u) = 6(u - x)(o) a + 9{x - u)(o) b + A§>6(x - u) + A$Sf(x - u) + . . . , (14) 

where 9{x) is the step function equal to 1 for x > and zero for x < 0, and the prime denotes 
differentiation. If p ab (u)du is the probability that the curve intersects the axis y = in the 
interval [u,u+du], with. p ab (u) = p a b{—u) and dup a b(u) = 1, then the average magnetization 
is 



(a(x, 0)) ab = / dup ab {u) a ab {x\u) + ... 

J — OO 

f'OO f'X 

= {a) a du p ab (u) + (a) b du p ab {u) + A^ b 'p ab (x) - A { Jp' ab (x) + . . . (15) 

J X J — oo 

where the dots in the first line stay for the contribution of multi-kink states. Comparison 
with (fTT|) and (fT3|) . as it is easily seen, shows correspondence between the two expansions and 
determines 

Pab(u) = e ~^/R f (16) 

A^=c^/m ab . (17) 

The last as well as additional terms in (115)) should be compared with those produced in ([6]) by 
further expansion around 9 = 9' = 0. The gaussian form of p ab (u) is model independent and 



4 At a later stage in the expansion thickness is generated also by multi-kink terms in the boundary state JTJ, 
which will produce a bundle of thin interfaces rather than just one (Fig. [2^). Multi-kink contributions to Z a b are 
suppressed at large R as e~ MR , M being the total mass. 
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of course coincides with the scaling limit of that deduced in [7] for the Ising models. As shown 
in the previous section the integral, non-local terms in (|15|) are entirely due to the pole term in 
(|12p . which in turn is produced by the non-locality of the kinks with respect to the spin fielcfl. 

For any lattice configuration there is a cluster (let us call it the left cluster) formed by 
nearest neighboring spins with value (color) a and whose external perimeter is formed by the 
left half of the boundary of the strip together with a path connecting the two boundary changing 
points. This path, whose identification may be ambiguous and require some lattice dependent 
prescription, becomes a simple curve in the continuum limit. A second curve connecting the 
two boundary changing points completes the perimeter of the cluster of color b anchored to the 
right half of the boundary (the right cluster). In general the two curves, which can touch but 
not intersect, enclose other clusters in between them. Among these intermediate clusters, those 
adjacent to the left (right) cluster have color different from a (b). The first few terms in (|14p and 
([15]) are compatible with a picture in which a uniform magnetization (a) a ((a)b) is assigned to 
the region enclosed by the perimeter of the left (right) cluster: the two curves become coincident 
as R/E, — > oo, with asymptotic passage probability density given by (fT6|) ; the first deviation from 
this situation as R/l; decreases is expected to happen via bifurcation and recombination around 
a cluster of color c 7^ a,b (Fig. [2b), and to be associated to the term containing . We will 
see in the next section that this term is indeed absent in the Ising model, where bifurcation is 
not allowecFI. 



4 Specific models 

Potts model. The lattice Hamiltonian [16] 

Hpotts = -J ^2 ^(xi),s(x 2 ). = l,...,q, (18) 

(xi,X2> 

is invariant under global permutations of the values of the spins s(x). For J > in two 
dimensions the phase transition is continuous for q < 4 and above J c there are q degenerate 
vacua located at the vertices of a hypertetrahedron in the (q-l)-dimensional order parameter 
space. Kinks with equal masses run along the edges of the hypertetrahedron and all the vacua 
are adjacent according to the definition given in section 2. The results we obtained for the 
magnetization profile apply to each component <r c (x) = 5 S ( X ) ]C — 1/q, c = 1, . . . ,q, of the spin 

J The gaussian probability density (|16p indicates an effective brownian behavior of the interface. The conver- 
gence of the interface to a brownian bridge for all subcritical temperatures has been proved in [TJ for the Ising 
model and in [15] for the g-state Potts model. We thank Y. Velenik for bringing these references to our attention. 

6 If we consider the form factor of the energy density e, the residue at iir is given by (O with a replaced by e, 
and vanishes because the expectation value of e is the same in all stable phases. This reflects the fact that even 
for a domain wall excitation the energy density is spatially localized (on the wall). 

7 Splitting into an odd number of paths (Fig. however, is allowed and encoded by subsequent terms in the 
expansion. 
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field; taking into account that Ylc=i °c = and {a a )b = qSab ^ l 1 (a a ) a , one obtains 



(cr c (x,0)) ab 



(<r a ) a 



q($ca + $cb) - 2 2q(5 ca - 5cb) 



dup(u) 



+ [2 - q(5 ca + 5 cb )] — p(x 
m 



? - 1 9-1 Jo 

+ . . . , (19) 



where p(u) is (|16p with m ab = m. Potts field theory is integrable [T7] and from the known form 
factors [12] we obtain 

B ~ 2^3' 3^3 ^ 



for q = 3, 4, respectively. For a,6 the integral term in (119]) is absent and the x-dependence 
of the magnetization profile is entirely due to the structure of the interface. The gaussian 
term in (|19p is produced by the leading deviation from the picture of the interface as a simple 
curve separating the phases a and b, which was argued in the previous section to correpond 
to the appearance of an island of the phase c via bifurcation and recombination of the curve. 
Bifurcation is not allowed in the Ising model, and indeed the coefficient of p(x) vanishes at q = 2, 
where c necessarily coincides with a or b. Directly at q = 2, the same conclusion is obtained from 
the explicit form F° ba {9) = i{a) a tanh((9/2) of the spin form factor, implying that (|12p contains 
only the terms with k odd. 

Percolation. It is well known [1H\ \TU\ that the partition function of the Potts model admits an 
expansion over Fortuin-Kasteleyn (FK) clusterqj of spins with the same color which as q — > 1 
become the clusters of random percolation. If we consider the Potts model with boundary 
conditions of type ab on the strip, the probability {6 s ( x ,o),a)ab that the spin s(x,0) has color a is 
given by the probability P(x, 0) that it belongs to a cluster touching the part of the boundary 
with x < (which has color a), plus 1/q times the probability 1 — P(x, 0) — P(—x,0) that it 
belongs to a bulk cluster. This can be rewritten as 

M*, 0))o6 = ?—^P(x, 0) - -P(-x, 0) . (21) 
q q 

The FK expansion of Z ab does not contain configurations with clusters connecting the boundary 
regions with x < and x > 0; this restriction is inherited by the percolation problem we 
consider. When the occupation probability p for the sites is above the percolation threshold 
p c (this corresponds to our case J > J c in the Potts model), even for R — > 00 there is a 
positive probability of having a spanning cluster which connects the upper and lower parts 
of the boundary with x < (Fig. [3]). Then the probability P(x,0)\ q= i that the site (x,0) is 
connected to the left part of the boundary is given by the probability P s (x,0) that it belongs 
to such spanning cluster plus the probability P ns (x, 0) that it belongs to a cluster touching only 
the upper or lower edge. Since the clusters of the latter type have an average linear extension 



The FK clusters differ from the 'geometrical' spin clusters we referred to in the previous section for the fact 
that nearest neighbors with the same color do not necessarily belong to the same cluster. 
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Figure 3: For percolation on the strip of width R, consider only the configurations without 
clusters connecting the part of the edges with x < to the part with x > 0. Eq. (I22p gives the 
probability that a point on the x axis belongs to a cluster s spanning between the negative part 
of the edges, in the scaling limit above p c and for R 3> £ = l/2m. 



of order 1/m, P ns (x,0) vanishes exponentially when mR — > oo. Hence in this limit we have 
P s (x,0) ~ P(x,0)\ q=1 , and from ([15]! . ([2I Jl we obtain 



P.(x,0) = - 

where P = lim g _^i -§j (cr a ) a is the probability that a site belongs to the infinite cluster on 
the infinite plane, and the mass m is related to the exponential correlation length above p c as 
£ = l/2m (see [191 120]); 7 = Res g= ii?(g) is not known because for q — > 1 the form factor © is 
available only for real rapidities |19| . 

Ashkin- Teller model. The lattice model corresponds to two Ising spins <ti(x), <72(x) = ±1 on 
each site interacting as specified by the Hamiltonian 

HAT = ~ ^2 { J [°'l( X l) 'l( X 2) + °"2(xi)cr2(x2)] + J4fTi(xi)cri(x 2 )cr2(xi)cr 2 (x2)} . (23) 

( x l, x 2> 

Each site can be in one of four states (01, 02) that we label a = 1, 2, 3, 4, corresponding to (+, +), 
(+,—), (—,—), (—,+), respectively. The model, that we consider for J > 0, is well known to 
possess a line of critical points parameterized by J4 |21| . In the scaling limit close to this line 
it renormalizes onto the sine-Gordon field theory (see [22\ I23j). where a parameter (3 2 plays the 
role of J4; /3 2 = 47r describes a free fermionic theory and corresponds to the decoupling point 
J4 = 0. Below critical temperature the model possesses four degenerate vacua |O a ) and for any 
value of /3 2 there are kinks |i^ a ,a±i(mod4)) with the same mass m. 

For f3 2 < 4tt the interaction among these kinks (which correspond to sine-Gordon solitons) 
is attractive and produces bound states, the lightest of which have mass m' = 2m sin 2{&f-0 i ) 
and are kinks \K a ^ a±2 (modA)) ■ Hence, in this regime any pair of vacua is connected by a single- 
kink excitation, all the vacua are adjacent and the boundary state \B a b) has in any case the 
form dU). The results of the previous sections apply with surface tensions ^ a ,a±i(mod4) = m 




7rmR 



e -2mx 2 /R _|_ 



(22) 
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and £ a ,a±2(mod4) = m '\ the bifurcation coefficients (fT7|) can be obtained from the form factors 
computed in [23]. For (3 2 = 2vr the masses m and m' coincide and one recovers the q = 4 
Potts model. 

For [3 2 > 4ir there are no bound states and the vacua with indices differing by two units are 
not adjacent, with surface tension £ a ,a±2(mod4) = 2m. In this case (pQ) is replaced by 



\B a , a±2 (x ;t)) = e- itH+ixoP 



(24) 



£ / ^^^±2(01^2)1^0(^1)^2(02)) + ... 

_c=a±l ^ 7F 

(indices are intended mod 4) and (|14p has to be replaced by a description in terms of two 
interfaces (Fig. [2b). 

Let us mention that studies of cluster boundaries at criticality can be found in particular in 
[24[ [25] for the Potts model and in |26[ [271 128] for the Ashkin- Teller model; we refer the reader 
to [5U] for results on cluster densities in critical percolation. 

In summary, in this paper we have shown how field theory naturally accounts for phase 
separation in two dimensions in the scaling limit below criticality. The derivation is simple but 
requires to take into account that the large distance properties are determined by kink (domain 
wall) excitations. In the last decades two-dimensional field theories with kinks have been studied 
within the non-lagrangian framework of integrable field theory based on asymptotic states, form 
factors and spectral sums; the derivation of section 2 exploits this framework but is general 
and does not require integrability. We have shown how it naturally leads to a gaussian passage 
probability density for the interface and accounts for the deviations from curvilinear behavior. 
The leading deviation is produced by branching of the inferface and has been determined exactly 
for the Potts model. Specializations of the formalism to percolation and to the Ashkin- Teller 
model have also been discussed. 
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